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ABSTRACT 

We present a power spectrum analysis of the final 2dF QSO Redshift Survey cata- 
logue containing 22652 QSOs. Utihsing the huge volume probed by the QSOs, we can 
accurately measure power out to scales of ~ 500/i~^Mpc and derive new constraints, 
at z ~ 1.4, on the matter and baryonic contents of the Universe. Importantly, these 
new cosmological constraints are derived at an intermediate epoch between the CMB 
observations at z ~ 1000, and local {z ^ 0) studies of large-scale structure; the av- 
erage QSO redshift corresponds to a look-back time of approximately two-thirds of 
the age of the Universe. We find that the amplitude of clustering of the QSOs at 
z ~ 1.4 is similar to that of present day galaxies. The power spectra of the QSOs at 
high and low redshift are compared and we find little evidence for any evolution in 
the amplitude. Assuming a A cosmology to derive the comoving distances, r(z), to 
the QSOs, the power spectrum derived can be well described by a model with shape 
parameter T — 0.13 ± 0.02. If an Einstein-de Sitter model r(z) is instead assumed, a 
slightly higher value of F = 0.16 ± 0.03 is obtained. A comparison with the Hubble 
Volume ACDM simulation shows very good agreement over the whole range of scales 
considered. A standard {flm = 1) CDM model, however, predicts a much higher value 
of r than is observed, and it is difficult to reconcile such a model with these data. 
We fit CDM model power spectra (assuming scale- invariant initial fluctuations), con- 
volved with the survey window function, and corrected for redshift space distortions, 
and find that models with baryon oscillations are slightly preferred, with the baryon 
fraction fltj/flm = 0.18 ± 0.10. The overall shape of the power spectrum provides a 
strong constraint on flmh (where h is the Hubble parameter), with flmh = 0.19±0.05. 

Key words: cosmology: observations, large-scale structure of Universe, quasars: 
general, surveys - quasars 



1 INTRODUCTION 

The large-scale structure of the Universe represents one of 
the most powerful discriminants between cosmological mod- 
els. QSOs are highly effective probes of the structure of the 
Universe over a wide range of scales and can trace clustering 
evolution over a look-back time which is 70-80 per cent of its 
present age. In the linear regime {W <,r <, lOOOh^^Mpc), 
they are clearly superior to galaxies as probes of large scale 
structure by virtue of both the large volumes they sample 
and their flat n{z) distribution. Accurately measuring the 



clustering over these scales bridges the gap between the clus- 
tering results on relatively small scales from galaxy redshift 
surveys out to scales previously only probed by microwave 
background anisotropy experiments. 

The power spectrum perhaps provides the most natural 
description of the matter fluctuations that comprise large- 
scale structure; for a Gaussian field, the amplitude of the 
Fourier modes provide a statistically complete description 
of the density perturbations. One potential problem of us- 
ing QSOs as probes of large scale structure is that they, like 
galaxies, are biased tracers of the mass density held. QSO 
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imaging experiments at redshifts out to z ~ 1 — 2 suggest 
that the QSOs reside in relatively modest (~ L*) galaxy 
hosts (Rix et al. 2001; Ridgway et al. 2001) and that the 
clustering environment of UVX QSOs are similar to opti- 
cally selected galaxies (Croom & Shanks 1999). This implies 
that at z ^ 1 — 2 QSOs and galEixies have a similar bias with 
respect to the underlying mass density field (Kauffmann & 
Haehnelt 2002). It is likely that the bias is scale-dependent, 
at least at small scales (Blanton et al. 1999) . At larger scales 
(r >> lOh^^Mpc), however, whore QSO clustering statis- 
tics prevail, the shape (although not the amplitude) of the 
measured power spectrum can reasonably be assumed to 
apply to the mass and galaxies as well, since it is unlikclj' 
that the QSO bias will be scale dependent in this regime 
(Coles 1993; Mann, Peacock & Heavens 1998). For the pur- 
poses of modelling the QSO P(k) in this paper we assume 
that any bias is scale independent on the scales considered 
(40 < r < SOOh-^Mpc). 

With the principal aim of determining the form of QSO 
clustering, and hence producing strong new constraints on 
cosmological models, we have used the 2dF multi-fibre spec- 
trograph on the Anglo- Australian Telescope to make a large 
QSO rcdshift survey The 2dF QSO Rcdshift Survey (2QZ) 
comprises two 5 deg x 75 deg declination strips, one at the 
South Galactic Pole and one in an equatorial region in the 
North Galactic Cap. QSOs are selected by ultra-violet ex- 
cess (UVX) in the u — bj : bj — r plane. The effectiveness 
of our QSO selection criteria has boon demonstrated with 
an overall QSO sky density of 35deg~^ to 6j < 20.85 being 
achieved, at a completeness of 93 per cent. Spectroscopic ob- 
servations arc now complete, and we have been able to iden- 
tify 86 per cent of the colour-selected candidates. Around 
half, 47 per cent, of the 47768 candidates are positively iden- 
tified as QSOs, leading to a final catalogue containing 22652 
QSOs. The QSOs typically have redshifts 0.3 <« < 3, and 
probe a volume of around 2 x 10® Mpc^ (assuming an 
Einstein-do Sitter cosmology). 

The spectra of the first 10000 2QZ QSOs were released 
in April 2001 (Croom et al 2001a) and can be obtained at 
http://www.2dfquasar.org. Hoylo et al. (2002) used these 
QSOs to measure the QSO power spectrum. They found that 
the QSO power spectrum has similar amplitude to that of 
local galaxies, and little evolution is seen with redshift. The 
main conclusion was that they detected large-scale power 
significantly in excess of the standard {flm = 1) CDM pre- 
diction, strongly ruling out this model and more consistent 
with CDM models with lower mass content and a non-zero 
cosmological constant (ACDM). 

In this paper, we ro-apply the power spectrum analysis 
described in Hoyle et al. (2002) to the final 2QZ catalogue 
containing 22652 QSOs. As well as having more than twice 
as many QSOs as available to Hoyle et al., this dataset has 
a much cleaner window function, significantly reducing the 
potential systematics on large scales. We briefly describe the 
2QZ data set in Section 2, including the angular and radial 
selection functions which are required to estimate the survey 
window function. Section 3 discusses the power spectrum 
estimation method, and the resulting QSO power spectrum 
is presented in Section 4. In this section we compare the QSO 
power spectrum to that of present day galaxies and clusters, 
and investigate evolution in the clustering amplitude. Then, 
in Section 5, we compare the QSO power spectrum with 



models of large scale structure, deriving strong constraints 
on the baryon and matter content of the Universe, before 
drawing conclusions in Section 6. 



2 THE 2QZ CATALOGUE 

The 2QZ comprises two 5 deg x 75 deg declination strips, one 
at the South Galactic Pole (centred at 5 = —30°, with 
2l''40™ ^ a <, 03''15™) and one in an equatorial region 
in the North Galactic Cap (centred at 5 = 0° with 09*' 50™ 
^ a <, 14*" 50™). We wiU refer to these regions as the SGC 
and NGC respectively. QSOs axe selected to bj < 20.85 by 
ultra-violet excess (UVX) in the u — bj : bj — r plane, at 
a completeness of 93 per cent. The spectra of the objects 
are obtained using the 2dF instrument on the AAT and re- 
duced using the 2dF pipeline reduction system (Bailey & 
Glazebrook 1999). Objects are identified as QSOs by an au- 
tomated procedure known as AUTOZ (Miller et al., in prep.). 
AUTOZ also determines the QSO redshifts; these have then 
been visually checked by two independent observers. The 
final 2QZ catalogue, containing 22652 QSOs (only those 
QSOs with quality 1 are used in this analysis; see Croom 
et al. (2001a) for details), covers an area of 740 deg^, and 
the data will be released to the community in the first half 
of 2003. 

2.1 The QSO selection function 

In order to estimate the QSO power spectrum, we need to 
take into account the various selection effects introduced 
when constructing the survey. This is achieved by generat- 
ing a catalogue of random points that mimics the angular 
and radial selection functions of the QSOs but otherwise is 
unclustered. 

The rapid luminosity evolution seen in QSOs (Boyle et 
al. 2000) means that at each redshift, over the redshift range 
we are considering, we are studying roughly the same part of 
the QSO luminosity function, relative to M* . In areas where 
the spectroscopic completeness is low, due for example to 
poor seeing conditions, we are sampling on average slightly 
brighter QSOs. This just has the effect of sampling QSOs 
further up the luminosity function at all redshifts, rather 
than sampling objects at lower redshift, as is typical with 
galaxy redshift surveys. Even if we cut back the magnitude 
limit of the survey by half a magnitude from 6j — 20.85 
to bj — 20.35, the mean redshift only drops slightly from 
z = 1.48 to 2 = 1.45. Therefore, for this work we have 
assumed that the radial and angular selection functions can 
be appUed independently, and consider each in turn below. 

2.1.1 The angular selection function 

The QSOs were observed simultaneously with galaxies from 
the 2dF Galaxy Redshift Survey (2dFGRS), and an adaptive 
tiling algorithm was developed to allow as many galaxies 
and QSOs as possible to be observed in each pointing. The 
pointings overlap in high density regions to maximise the 
coverage in all areas of the survey, and with this algorithm, 
we have achieved a completeness of 93 per cent. In addition 
to the 7 per cent of objects on which we were unable to 
place fibres, regions around bright stars were omitted from 
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the input catalogue, and hence the angular mask of the 2QZ 
is quite complicated. 

The survey also suffers from spectroscopic incomplete- 
ness. We have failed to confirm the identity of approximately 
14 per cent of the objects observed. The spectra of these 
objects are typically very noisy, with a lack of discernible 
spectral features. As some classes of objects are easier to 
identify in low S/N spectra than others, due to, for exam- 
ple, distinctive emission lines, the spectroscopic complete- 
ness of QSOs in the survey may differ from that of stars, or 
the input catalogue as a whole. To correct for spectroscopic 
incompleteness we therefore need to estimate the fraction 
of QSOs that remain in those unidentified spectra. By ex- 
amining the results from objects that have been observed 
more than once, such an estimate can bo obtained. For ob- 
jects unidentified after their first observation, but identified 
in subsequent observations, the relative fractions of QSOs 
and stars are remarkably similar to those of the identified 
objects in the full survey. Hence we estimate that QSOs 
account for approximately the same proportion (~ 59 per 
cent) of the identified and unidentified populations and so 
the spectroscopic incompleteness of QSOs is in fact the same 
as that of the input catalogue as a whole. However, the rate 
of identifications obtained through the re-observation of pre- 
viously unidentified objects is significantly lower than the 
survey average, and hence it is still possible that the frac- 
tion of QSOs remaining in the unidentified population could 
bo somewhat lower. Through consideration of the fraction 
of QSOs unidentified after their first observation that are 
subsequently identified, we can set a lower limit to the frac- 
tion of QSOs in the remaining unidentified population at 
about 40 per cent. We will consider the effect of different 
spectroscopic completeness corrections in Section 4.1. 

To match the angular selection function of the QSOs, 
we construct a completeness map, estimating completeness 
both to observational and spectroscopic effects in regions 
defined by the intersection of 2dF pointings. In each such 
region we calculate the fraction of objects contained in the 
input catalogue that have been identified. As the spectro- 
scopic incompleteness of QSOs is approximately the same 
as that of the full input catalogue, the proportion of QSOs 
in unobserved objects should be the same as the proportion 
of QSOs in observed, but unidentified objects. Hence we can 
take this fraction to be the overall (observational and spec- 
troscopic) completeness in that region. As we are applying a 
spectroscopic completeness correction, wc have decided not 
to apply simultaneously a minimum spectroscopic complete- 
ness cut, as this would reduce the size of our sample. 

We also account for the extinction duo to Galactic dust, 
using the estimate of the dust reddening, E{B — V), as a, 
function of position given by Sclilcgcl, Finkbeiner & Davis 
(1998). As Galactic dust changes the effective magnitude 
limit, we weight the random distribution by 

M/ext(a,<5) = 10-''Vi (1) 

where Abj = 4:.035E{B - V) and /3 = 0.3, the slope of the 
integral QSO number-magnitude relation at the magnitude 
limit of the survey, f)j=20.85. As we apply a dust correction 
to the angular mask, we need to correct for dust extinction 
in the observed radial selection function (below) to avoid 
including the effect of dust twice. 
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Figure 1. The radial number density of QSOs in bins of Ar = 
SOfe^^Mpc in tlie NGC (solid line) and SGC (dashed hnc), cor- 
rected for tlie average dust extinction in each region and calcu- 
lated assuming an EdS r{z) (top) and the A r{z) (bottom). Over- 
laid are the 10th order polynomial fits used to create the random 
catalogues, plotted over the redshift range 0.3 < z < 2.2. The 
number density varies very slowly as a function of scale; almost 
constant, assuming the EdS r{z) and decreasing only by a factor 
of approximately three over the range 0.3 < 2 < 2.2, assuming the 
A r{z). At higher redshifts, the number density decreases rapidly. 

2.1.2 The radial selection function 

We arc measuring QSO clustering as a function of comov- 
ing distance, and so we need to adopt a cosmology to con- 
vert from redshift into comoving distance. In this paper we 
consider an Einstein-de Sitter {Qrn~1.0, f2A=0.0) cosmol- 
ogy r{z) (EdS hereafter) and an ilm^O.S, Oa=0.7 cosmology 
r{z) (A hereafter). 

Figure 1 shows the radial selection function of QSOs in 
the NGC (solid line) and SGC (dashed line), corrected for 
the average dust extinction in each region and calculated 
assuming the EdS r{z) or the A r{z). Due to the contri- 
bution from the host galaxy, QSOs at low redshift are lost 
from our input catalogue because they do not appear stellar. 
The host galaxy also reddens faint QSOs causing then to be 
preferentially lost from the sample. To limit this incomplete- 
ness, we restrict our analysis to z > 0.3 (r > 740/i^^Mpc, 
assuming EdS, or r > 835/i~^Mpc, assuming A). At high 
redshifts the UVX colour selection technique breaks down 
as the Lyman-alpha forest enters the U band, and the com- 
pleteness of the survey rapidly drops. Therefore, wo also 
limit our analysis to 2; < 2.2 (r < 2645/i~^Mpc, assuming 
EdS, or r < 3820/i~^Mpc, assuming A). Over the redshift 
range 0.3 < z < 2.2 the completeness of the 2QZ colour 
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Table 1. The radial selection function is fitted separately in 
NGC and SGC, assuming either the EdS r{z) or the A r{z) to 
derive comoving distances from the QSO redshifts. This table 
shows the values of the 10 coefficients obtained from a 9th or- 
der polynomial fit, N{r) = ^2i—o ( lOOO *° ^^'^ number of 
QSOs in each strip, in the range 500 < r < 4200/i^^Mpc (A), or 
500 < r < 2800/i-lMpc (EdS), in bins of Ar = 50/i-lMpc. The 
QSO number density, n(r), is then obtained via n(r) = ■^N{r), 
where aj^^^ = 0.2456, and OgGC ~ 0.2281 are constants that 
take into account the average completeness and extinction in each 
strip. 





NGC, EdS 


NGC, A 


SGC, EdS 


SGC, A 


Co 


104.623 


47.1019 


295.615 


124.692 


Cl 


-448.886 


-216.212 


-987.737 


-440.183 


C2 


576.207 


339.671 


615.208 


463.583 


C3 


-35.6319 


-154.792 


883.628 


-44.2143 


C4 


-176.484 


18.5957 


-858.731 


-92.7632 


C5 


17.0201 


0.628647 


-20.3508 


20.7977 


Ce 


48.8416 


1.83275 


188.357 


8.12686 


C7 


-13.6702 


-0.301663 


-13.6225 


-2.35053 


Cs 


-0.639615 


-0.126688 


-21.5386 


-0.123221 


Cg 


0.269846 


0.022045 


4.13687 


0.049514 



selection technique is J; 90 per cent. The SGC has a slightly 
higher normalisation (15 per cent) than the NGC, even after 
correcting for the higher levels of dust in the NGC, but very 
similar shape, except at low redshift (r < 1200/i~^Mpc, or 
z < 0.45) where the SGC suffers from higher incompleteness, 
probably due to star/galaxy separation differences. For this 
reason we treat the radial selection function separately in 
the two strips. The QSO number density varies very slowly 
as a function of redshift. Assuming EdS it is almost con- 
stant over the range 0.3 < z < 2.2, whereas assuming A, it 
decreases only by a factor of approximately three over the 
range 0.3 < 2 < 2.2. Applying these redshift cuts therefore 
allows us to consider clustering statistics without the need 
for a weighting scheme. This gives us a sample of 8704 QSOs 
in the NGC, and 10845 QSOs in the SGC. 

We estimate the space density of QSOs as a function 
of comoving distance, n(r), via a 10 coefficient (9th order) 
polynomial fit to the number distribution of QSOs, N{r), 
in each strip in bins of Ar = 50ft~^Mpc. Further details 
of the fits, including the coefficients, are given in Table 1. 
The resulting n(r) distributions, used to create the ran- 
dom catalogue, are plotted in Fig. 1 over the redshift range 
0.3 < z < 2.2. A low order polynomial would provide a 
poor fit to the n(r) distribution. As the order is increased, 
the quality of the fit, as measured by the statistic, im- 
proves dramatically until approximately an 8/9 coefficient 
polynomial is used. Further increasing the order of the fit, 
however, hardly reduces the value. If a much higher order 
polynomial wore used, it may also start to fit real features in 
the n(r) distribution, and hence artificially remove power. 
Hence a 10 coefficient polynomial is chosen so as to provide 
an adequate fit to the n(r) distribution without over-fitting. 

A poor fit to the n(r) distribution would artificially 
introduce waves in the normalised QSO distribution, and 
hence add power on the largest scales. Indeed, when the A 
power spectrum is measured using a radial selection function 
derived from a 5 coefficient polynomial fit, and comparing to 



the result using a 10 coefficient polynomial fit, we find that 
the power is significantly enhanced on scales > 200/i^^Mpc 
(log(A;//iMpc^^) < —1.5). Slightly increasing the order of 
the polynomial fit, however, has little effect our power spec- 
trum estimate. Rcmeasuring the A power spectrum using a 
15 coefficient polynomial fit, and comparing to the result us- 
ing a 10 coefficient polynomial fit, we find that out to scales 
of 400/i~^Mpc (Iog(fc/;iMpc"^) > -1.8), the power spectra 
agree in amplitude to within 4 per cent, considerably smaller 
than the fractional error at this scale. 

2.2 Mock 2QZ catalogues 

To test our method of power spectrum estimation described 
below, as well as our error estimates, we have used the 
huge Hubble Volume ACDM simulation (Fronk ot al. 2000, 
Evrard et al. 2002), where the particles have been output 
along an observer's past light cone to simulate the 2QZ. 
The cosmological parameters of the simulation are fib =0.04, 
ficDM=0.26, fiA=0.7, Ho=70 km s"^Mpc"^ and ag = 0.9. 
One billion mass particles arc contained within a cube that 
is 3,000ft^^Mpc on a side. To create realistic QSO mock 
catalogues to match the final 2QZ sample, we bias the mass 
particles to give a similar clustering pattern, and introduce 
the same angular and radial selection function as seen in 
the final 2QZ sample, described in the previous section. For 
further details see Hoyle (2000), where the Hubble Volume 
simulation was compared to the 10k QSO catalogue. 



3 POWER SPECTRUM ESTIMATION 

The power spectrum estimation is carried out as described 
in Hoyle et al. (2002), using the method outlined in Tadros 
& Efstathiou (1996). The two declination strips are treated 
separately. They arc embedded into a larger cubical vol- 
ume, and the density field is binned onto a 256'* mesh, using 
nearest grid-point assignment. The power spectrum of each 
region is estimated using a Fast Fourier Transform (EFT), 
and the average of the resulting power spectra is taken. The 
Fourier modes arc binned up logarithmically to reduce the 
covariance between each bin. The EFT is unreliable at small 
scales; for a 256^ EFT and a box of size 4000/i"^Mpc, the 
Nyquist frequency of the transform is /cNyq = 0.2/i,Mpc~^. 
The limit out to which the power spectrum measurement 
is reliable depends on the grid assignment scheme, and is 
a fraction of fcNyq- Using the nearest grid point assignment 
scheme, fcum « 1/2 AiNyq (Hatton 1999), which, for the above 
box size, corresponds to a value of log(fciim//i Mpc"'^) w — 1 
or a scale of ~ 60/i^^Mpc. 

The measured power spectrum is convolved with the 
power spectrum of the window function. As the QSO co- 
moving number density is almost constant out to very large 
scales, we arc effectively considering a volume limited sam- 
ple. Hence each QSO carries equal weight, and the survey 
window function, W{x.), simply takes a value of unity in the 
volume of the universe included in the survey, and zero else- 
where. This is approximated using a catalogue containing a 
large number of unclustcrcd points with the same radial and 
angular distribution as the survey, and its power spectrum 
is calculated in a similar manner to that of the data. 

To enable measurements of the power spectrum on 
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smaller scales, we use a direct method of computing the 
Fourier transform. This method is time intensive and can 
only be applied to a small number of particles (so cannot be 
applied to the random catalogue). Therefore, we make the 
assumption that the window function has a small amplitude, 
and so negligible effect on the power spectrum estimate. Us- 
ing mock catalogues drawn from the Hubble Volume simular- 
tion, Hoyle (2000) demonstrated that on scales smaller than 
w 100/i"^Mpc (log(fc/ftMpc~^) ?a -1.2) this approximation 
is accurate. 

Finally, to obtain a single 2QZ power spectrum for each 
assumed cosmology, we average the power spectrum of the 

QSOs in the NGC and SGC strips together, weighting by 
the inverse of the variance on each scale. 



3.1 The window function 

In Fig. 2 we show the window function, IWfep, of the 2QZ, 
measured assuming either EdS or A. The window function 
is anisotropic, due to the shape of the survey and varia- 
tions in completeness, described by the selection function. 
In order to compare model power spectra with that deter- 
mined here for the 2QZ, we need to convolve the model 
power spectra with the window function. In the case where 
the power is isotropic, the same result would be obtained by 
convolving with the spherically averaged window function. 
In reality, whilst we expect the real-space power spectrum 
to be isotropic, redshift-space distortions introduce small 
anisotropics into the shape of the redshift-space 2QZ power 
spectrum (see Outram et al. 2001). However, the Hubble Vol- 
ume simulations suggest that these anisotropies have little 
effect on the spherically averaged power spectrum estima- 
tion, and so we only consider the spherically averaged win- 
dow function from here on. We also show an analytic fit to 
the spherically averaged window function in Fig. 2: 

< \Wkf >Eds= [1 + (fc/0.00204)^ + (fc/0.00362)*]-^ (2) 

< \Wkf >A= [1 + (fc/0.00132)^ + (fc/0.00258)^]"' (3) 

Due to the huge volume probed by 2QZ, the window 
function is much more compact than that of galaxy surveys, 
and it takes the form of a steep power law, proportional to 
fe"* out to scales of w 600/i"^Mpc (log(fe//iMpc"^) w -2) 
for A. For comparison we also show an analytic fit to the 
2dFGRS window function (Percival et al. 2001) in Fig. 2. 
The 2QZ window function has a very similar shape to that 
of the 2dFGRS, but shifted to larger scales; a factor of 2.5 x 
larger for the 2QZ EdS window function, or a factor of 3.5 x 
for A. 

Fig. 3 demonstrates the effect that the 2QZ window 
function has on power spectra. On scales < 200/i~^Mpc 
(log(fe//iMpc~^) > —1.5), the window function has very 
little effect on the true power spectrum, whilst on scales 
larger than this the true power is slightly suppressed. Even 
out to scales of fa eOOfe^^Mpc (log(fe//i Mpc"^) ^ -2), the 
effect of the window function is considerably smaller than 
the statistical errors on the power spectrum estimate. The 
broader 2dFGRS window function heis a much larger effect, 
not only suppressing power on scales > lOO/i^^Mpc, but 
also smoothing out any baryonic 'wiggles' at smaller scales 
(Miller, Nichol & Chen 2002), and introducing covariance 
between the P{k) data points. 
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Figure 2. The window function of the 2QZ NGC region, mea- 
sured assuming EdS (top) and A (bottom). The window function 
is estimated using a random catalogue created to match the 2QZ 
selection function, containing 25 times as many points as the QSO 
catalogue. The points show measurements in individual k modes 

from the FFT grid. The tail to high k at log(|W^fcp) 5 is due 

to shot noise in the window function estimate. The bold points 
show the binned, spherically averaged window function, and the 
solid line is an analytic fit to the spherically averaged window 
function. In either cosmology, the window function for the SGC 
is almost identical to that of the NGC. Due to the huge volume 
probed by 2QZ, the window function is much more compact than 
that of galaxy surveys, and has the form of a steep power law 
out to scales of 600/1"^^ Mpc. For comparison an analytic fit to 
the spherically averaged 2dFGRS window function is also shown 
(dashed line) , together with the same window function shifted to 
larger scales by a factor 2.5 x (EdS) or 3.5 x (A) (dotted line). 



3.2 Error determination 

The errors on the power spectrum determination are es- 
timated using the method of Feldman, Kaiser & Peacock 
(1994; FKP), equation 2.3.2. The error is given by 



P2(fc) 



VkVs 



(4) 



where Vk is the volume of each bin in fc-space, estimated by 
Vk = ATfc ( Afc)^ with Nk the number of independent modes 
in the A;-shell and (Afe)^ the volume of one fe-mode. Vs is 
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Figure 3. The effect of convolving a model power spectrum with 
the 2QZ A window function. The solid line shows a linear ACDM 
power specrum, with Clm = 0.3, Hq = 70kms~^Mpc~'^, and Clf, = 
0.04. The dotted line shows the same power spectrum convolved 
with the 2QZ A window function, and the dashed line shows the 
power spectrum convolved with the 2dFGRS window function. 



the volume of the survey. (We have assumed that the QSOs 
carry equal weight when estimating the power spectrum and 
associated errors). Hoyle (2000) demonstrated that over a 
wide range of scales these errors were found to be in rea^ 
sonable agreement with errors estimated from the disper- 
sion over power spectrum estimates from the Hubble Volume 
mock catalogues. On the smallest scales, non-linearities may 
lead to a larger error than the FKP error estimate (Meiksin 
& White 1999). 



4 THE QSO POWER SPECTRUM 

In Figure 4 we show the power spectrum estimate obtained, 
using the method outlined above, for the 2QZ dataset re- 
stricted to QSOs with redshift in the range 0.3 < z < 2.2. 
We compare the NGC and SGC power spectra measured 
assuming the EdS (a) or A (b). The combined power spec- 
tra are also shown, assuming EdS (c) or A (d). The power 
spectra measured from the two strips agree well over a wide 
range of scales, with a dispersion that is consistent with our 
error estimate. The errors shown here and throughout the 
paper are la errors estimated using the FKP method. The 
power spectrum estimates obtained are also given in Table 2 
(assuming EdS) and Table 3 (assuming A). 

We compare this power spectrum estimate (assum- 
ing A) to that obtained by Hoyle et al. (2002) using 
the 10k 2QZ catalogue in Fig. 5. On scales smaller than 
log(A://i Mpc^^) ~ —1.6, the spectra are very similar, with 
the main difference being the reduced errors for the new esti- 
mate due to the increased sample size. On large scales, how- 
ever, the 10k power spectrum estimate is slightly higher than 
the latest determination. This over-estimate was largely due 
to problems accounting for the complicated window of the 
incomplete 10k survey. However, with this difference in mind 
we now examine some of the other potential contributions 



Table 2. Values of the QSO power spectrum P{k) estimated 
from QSOs with redshift in the range 0.3 < z < 2.2 contained 
in the 2dF QSO Redshift Survey Catalogue, assuming an EdS 
cosmology r{z). The errors are estimated using the method of 
Feldman, Kaiser & Peacock (1994). 



log(fe/hMpc-i) 


P{k)/h-^ Mpc3 


AP{k)/h-^Mpc^ 


-1.91250 


7497.86 


15528.9 


-1.83750 


9975.61 


14241.7 


-1.76250 


11055.6 


9826.21 


-1.68750 


14009.9 


8392.96 


-1.61250 


17765.4 


6605.29 


-1.53750 


14482.3 


4921.19 


1 AfiOr,(\ 
- J. .^uzou 


77r:Q on 
/ ( oo.oy 


Qr:i e 1 n 

OO-LO. -LU 


-1.38750 


10000.78 


2758.25 


-1.31250 


7918.47 


2079.62 


-1.23750 


10850.0 


1668.21 


-1.16250 


8722.57 


1252.32 


-1.08750 


6644.70 


940.720 


-1.01250 


4991.80 


710.215 


-0.937500 


3716.06 


538.652 


-0.862500 


2851.89 


411.061 


-0.787500 


1703.66 


312.109 


-0.712500 


1924.72 


241.619 


-0.637500 


1194.02 


184.568 


-0.562500 


900.635 


141.926 



Table 3. Values of the QSO power spectrum P{k) estimated 
from QSOs with redshift in the range 0.3 < z < 2.2 contained in 
the 2dF QSO Redshift Survey Catalogue, assuming a A cosmol- 
ogy r{z). The errors are estimated using the method of Feldman, 

Kaiser Peacock (1991). 



log(fe/hMpc-i) 


P(fe)//i-3 Mpc3 


AP(fe)//i-3 Mpc3 


-1.98750 


34308.1 


40765.3 


-1.91250 


31228.7 


30009.9 


-1.83750 


30285.0 


22002.5 


-1.76250 


66904.5 


20668.8 


-1.68750 


42073.0 


14148.3 


-1.61250 


27967.6 


10203.46 


-1.53750 


36571.7 


8169.82 


-1.46250 


19371.4 


5884.91 


-1.38750 


22035.7 


4602.21 


-1.31250 


20998.8 


3522.13 


-1.23750 


16388.0 


2646.94 


-1.16250 


16152.2 


2040.24 


-1.08750 


10279.27 


1531.78 


-1.01250 


8100.51 


1171.00 


-0.937500 


5352.25 


891.509 


-0.862500 


5013.59 


687.006 


-0.787500 


3816.35 


527.254 


-0.712500 


2233.81 


403.834 


-0.637500 


841.037 


309.547 


-0.562500 


800.378 


238.849 



to systematic errors in the power spectrum determination 
on large scales. 
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Figure 4. In Figures (a) and (b) wc compare the power spectrum from the NGC (fiUcd circles) and SGC (open 
circles). Both sets of points have been offset by 0.01 in log(fe/feMpc^-'-) for clarity. In Figures (c) and (d) we show 
the combined QSO power spectra. The EdS cosmology r{z) is assumed for (a) and (c) and A for (b) and (d). 



4.1 Large-scale systematic errors 



One possible source of systematic errors is the uncertainty 

in the zero-point cahbration of the UKST field plates that 
were used to define the input catalogue for the 2QZ sur- 
vey. The fields are 5degx5deg, and so any differences in 
the plate calibration could introduce power on the largest 
scales probed. The recent release of the SDSS EDR data 
(Stoughton et al. 2002) has given us the opportunity to ap- 
ply a new calibration to our data in the NGP to tost this 
effect, after converting the Sloan colours to our bands (Blan- 
ton et al. 2001; Fukugita, Shimasaku, & Ichikawa 1995). A 
comparison of the two calibrations suggests that our cal- 
ibration introduces platc-to-platc uncertainty in the QSO 
number density of the order of 4 per cent. In Fig. 6 we show 
the effect of using the SDSS calibration in the NGP on our 
power spectrum estimate. The difference in logP(A;) deter- 
mined using the two calibrations is shown by the solid line. 
On scales log(fc//t Mpc^^) > —1.8 there is virtually no dif- 
ference. On large scales, using the SDSS calibration would 
lead to a small decrease in the measured power spectrum. 
Unfortunately we have no similar data to re-calibrate the 



SGP strip, and so, because this effect is very small, we have 
continued to use our original calibration. 

A second potential source of systematic error is Galac- 
tic dust. We account for dust extinction using the estimates 
given by Schlegel, Finkbeiner & Davis (1998). If dust extinc- 
tion is not correctly accounted for then this could introduce 
large scale variations in the observed QSO number density, 
and hence power. The short-dashed line in Fig. 6 shows 
a comparison between the power spectrum estimates with 
and without a dust correction. Applying the correction sig- 
nificantly reduces power on scales log(fc//i Mpc^^) < —1.8. 
Whilst this gives us confidence in the dust maps, their reso- 
lution is fairly low, and so dust could still affect our estimate 
on smaller scales. 

Finally wc consider the effect of different choices of spec- 
troscopic incompleteness corrections. In the 10k power spec- 
trum analysis, only regions with spectroscopic completeness 
in excess of 85 per cent were considered, and no further cor- 
rection was applied to account for any QSOs still missing in 
the unidentified spectra. As identification is dependent on 
the quality of the spectra, which varies with each pointing 
due to observing conditions, the fraction of missing QSOs 
could vary slightly from region to region even in this 85 
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Figure 5. A comparison between tlie 2QZ power spectrum (as- 
suming A) from the final sample (filled circles) with that derived 
from the initial 10k catalogue (open circles; Hoyle et al. 2002). 
Both sets of points have been offset by 0.01 in log(A;//iMpc~^) 
for clarity. There is very good agreement, except on the largest 
scales where the excess power was introduced by the complicated 
incomplete window function in the 10k sample. 




log(k/h Mpc-') 



Figure 6. The effect of meiking several changes to the observa- 
tion mask are demonstrated. The solid line shows the change in 

the estimated power spectrum when SDSS photometry is used 
to calibrate the UKST plates in the NGP, rather than our own 
calibration. The short dashed line shows the effect of not correct- 
ing for dust. The long dashed line shows the difference between 
the power spectrum determinations using two different spectro- 
scopic completeness estimates. The dotted line shows the effect of 
including a spectroscopic completeness estimate when consider- 
ing only a sample with spectroscopic completeness > 85 per cent. 
For comparison, the double lines show the Icr errors for the power 
spectrum estimation. 

per cent complete sample, introducing large scale power. To 
investigate this we compare the power spectra determined 
using an 85 per cent spectroscopic completeness cut, then a 
further correction for observational completeness (as carried 
out in the 10k analysis) to one with the same 85 per cent 
cut, but also correcting for spectroscopic completeness, us- 



ing the same form of correction as in this paper; assuming 
that fraction of QSOs in unidentified objects is the same as 
the fraction in the identified spectra. The result is shown 
by the dotted line in Fig. 6. The lack of a spectroscopic 
completeness correction introduces a slight over-estimate in 
the power, especially at large scales. This may also partly 
explain the difference between the 10k and final power spec- 
trum estimates. 

The spectroscopic completeness correction we apply 
simply assumes that the spectroscopic incompleteness of the 
QSOs is the same as that of the full input catalogue, and 
hence that ~ 59 per cent of unidentified objects are QSOs. 
Whilst this is favoured by evidence from rc-obscrvations, 
there is still uncertainty in the proportion of QSOs remain- 
ing in the unidentified spectra, and a slightly lower fraction 
of unidentified QSOs is still compatible with the observa- 
tions. To test the effect of this uncertainty in the spectro- 
scopic completeness correction wc can construct a now com- 
pleteness map where wo instead assume that a lower fraction 
of approximately 40 per cent of objects that remain uniden- 
tified after observation are QSOs. A comparison between 
the power estimated using the two completeness corrections 
is shown by the long-dashed line in Fig. 6. The difference 
in power is very small on all scales, indicating that the un- 
certainty in the unidentified QSO fraction has a negligible 
effect on the power spectrum estimate. 

It is also possible that there are spectroscopic complete- 
ness variations across individual 2dF fields, which could in 
principle introduce power on smaller scales. By artificially 
introducing completeness gradients across the 2dF fields in 
the random catalogue, we can investigate the likely size of 
any systematic uncertainty due to this effect. Even with a 
completeness gradient of 20 per cent from the centre to the 
edge of each field, the change in the measured power spec- 
trum is no larger than that observed for any of the other 
potential systematics discussed above. As this artificially 
introduced gradient is much larger than the expected size 
of completeness variations across 2dF fields, given that the 
overall spectroscopic completeness for the survey is 86 per 
cent (and this is largely due to other effects such as weather 
conditions and seeing), we conclude that the power added 
by such an effect would be small, and we do not attempt to 
correct for this. 

All of the potential systematics discussed above are 
much smaller than the statistical errors on all scales. With 
the exception of the dust correction, they are of the order of 
AlogP(fc) ~ 0.05 even on scales of log(fc//iMpc"^) ~ -2. 
Therefore we have confidence that our estimate is robust out 
to these scales. 

4.2 The shape of the QSO power spectrum 

In CDM models, the power spectrum should turn over at 
a scale which depends on the shape parameter, F, of the 
power spectrum. The shape of the power spectrum and the 
position of the turnover contains information on cosmologi- 
cal parameters that can lead to strong constraints. In pure 
CDM models (i.e. without baryons) with a scale-invariant 
spectrum of initial fluctuations F = Q,mh. However, the in- 
clusion of baryons, or a tilted spectrum of initial fluctua- 
tions complicates the picture, so a measurement of F alone 
cannot be inverted reliably to give Qmh (Sugiyama 1995, 
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Figure 7. Model power spectra with varying shape parameter, 
r, are compared to the QSO power spectra for the two choices 
of cosmology. In the upper plot, the solid points show the QSO 
power spectra with the EdS ri^z). Overlaid is the best-fitting V 
model, with F = 0.16 (solid line). In the lower plot, the solid 
points show the QSO power spectra assuming the A Overlaid 
is the best-fitting V model, with F = 0.13 (solid line). On both 
plots, the dotted line shows a F = 0.1 model, the short dashed 
line shows a F = 0.2 model, and the long dashed line shows a 
F = 0.5 model. The latter model is favoured by a Standard CDM 
model, but clearly rejected by the data. 

Peacock 2000). For example, the A model we have assumed 
(f2m=0.3, f2A=0.7) to calculate the r{z) would have a value 
of r = 0.21, if wo further assume that ^ 70kms~^Mpc~^ 
(e.g. Tanvir et al. 1995; Mould et al. 2000), in the absence of 
baryons. However, when baryons are taken into account this 
model would then have an effective value of F = 0.17. For a 
full comparison, we also need to consider the effects of the 
window function, and rcdshift-spaco distortions. Wc have 
had to estimate the real-space distances to the QSOs from 
their observed redshift, and peculiar velocities, large-scale 
flows, geometrical assumptions and uncertainties in redshift 
measurements can all affect the shape and amplitude of the 
power spectrum wo have determined. Despite this, F is a 
very useful parameter for describing the general shape of 
power spectra, and allowing a simple comparison between 
different measurements or models. 



To parametrize the QSO power spectra we have mea- 
sured, we fit model CDM power spectra with different values 
of the shape parameter F. We assume the fit to the CDM 
transfer function given by Bardeen et al. (1986). For each 
value of F, we allow the amplitude to vary freely, and mini- 
mize the of the fit. The best fit models are shown in Fig. 7. 
For an EdS r{z), we find that models with F = 0.16 ± 0.03 
(IfT errors) provide the best fit to the power spectrtim, 
whereas assuming a A r{z), the power spectrum is slightly 
steeper, with F = 0.13 ±0.02. These low measurements of F 
indicate significant large-scale power, considerably in excess 
of Standard (flm = 1) CDM model predictions, and consis- 
tent with the A model prediction discussed above. In the 
case of the A r{z), whilst the power spectrum flattens at the 
largest scales probed, we see no evidence of the expected 
turnover below scales of ~400/i~^Mpc (log A: ~ —1.8). How- 
ever, in the slightly flatter EdS r{z) power spectrum there 
is a hint of a possible turnover by scales of ~400ft^^Mpc. In 
either case we can strongly rule out a turnover in the power 
spectrum at scales below ~300/i~^Mpc (log A: ~ —1.7). We 
will consider the effect of baxyons, the window function, and 
redshift-space distortions in Section 5 in order to derive con- 
straints on fim/i from the shape of the power spectrum. 

4.3 Comparison with P(k) from galaxy and 
cluster surveys 

We compare the QSO power spectrum with the power spec- 
tra of local galaxies, determined from the 2dFGRS (Pcrcival 
et al. 2001), PSCz (Sutherland et al. 1999), APM (Gaz- 
tanaga & Baugh 1998), and DUKST (Hoyle et al. 1999) 
galaxy surveys and the power specta of local Abell/ACO 
clusters (Miller & Batuski 2001) and APM clusters (Tadros, 
Efstathiou & Dalton 1998) in Fig. 8. Assuming a A r{z), 
the amplitude of QSO clustering at a 1.4 is similar to 
that of local galaxies, and almost an order of magnitude 
below that of galaxy clusters. Assuming an EdS r{z), the 
amplitude of QSO clustering is slightly lower than the local 
galaxy clustering amplitude. The A r{z) QSO power spec- 
trum appears steeper than that of the 2dFGRS, with more 
power on large scales. This is largely due, however, to the 
effects of the window function on the 2dFGRS power spec- 
trum, as demonstrated in Fig. 3. For a true comparison of 
the relative shapes of the galaxy and QSO power spectra, 
we will fit models convolved with their respective window 
functions in Section 5. 



4.4 The evolution of QSO clustering 

To investigate the effects of evolution on the amplitude of 
QSO clustering we split the QSO sample into two, a high 
redshift sample with z > 1.4, and a low redshift sample with 
z < 1.4, and compare their respective power spectra. The 
results are shown in Fig. 9. Assuming either EdS or A, we see 
virtually no evolution in the amplitude of clustering between 
the two samples, at 0.3 < 2; < 1.4 and lA < z < 2.2. 
To quantify this, we fit the amplitude of the best fitting 
model power spectra for each cosmology found in section 4.2 
(F = 0.16 for EdS and F = 0.13 for A) separately for each 
redshift bin. The results are shown in Table 4. We find that 
the amplitude of QSO clustering is almost identical in the 
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Figure 8. A comparison of galaxy and QSO power spectra. The 
EdS i'{z) is assumed in the upper plot and the A r(z) is assumed in 
the lower plot. The solid line shows the QSO power spectra. The 
errors are la FKP errors. The points show the power spectra of 
local galaxies, determined from the 2dFGRS (Percival et al. 2001), 
PSCz (Sutherland et al. 1999), APM (Gaztanaga & Baugh 1998), 
and DUKST (Hoyle et al. 1999) galaxy surveys and the power 
spectra of local Abell/ACO clusters (Miller & Batuski 2001) and 
APM clusters (Tadros, Efstathiou & Dalton 1998) 



two rcdshift bins. Assuming EdS, the difference in amplitude 
between the two redshifts is insignificant. Assuming A, we 
marginally detect a slight fall in amplitude of A log P = 0.05 
from high redshift to low redshift, at about la significance. 
This is in very good agreement with the results of the 2QZ 
correlation function analysis of Groom et al. (2001b). 

We fixed the value of F for the above analysis, as the 
shape parameter is assumed not to evolve with redshift (in 
the linear regime), in order to investigate the evolution of the 
amplitude of QSO clustering. We can test this assumption 
by comparing the value of F in the low and high redshift 
bins. For the A r{z), we find that this is indeed the case, 
with best fit values of F = 0.14 ± 0.03 at 0.3 < z < 1.4 
and r = 0.13 ± 0.03 at lA < z < 2.2. For the EdS r{z), 
however, there is marginal evidence for a slight flattening of 
the power spectrum at higher redshift, with best fit values 
of r = 0.15 ± 0.05 at 0.3 <« < 1.4 and F = 0.20 ± 0.05 



Table 4. The amplitude of model power spectra fitted to QSO 
power spectra in redshift bins, assuming either EdS or A. In both 
cases, the shape parameter, F, was fixed at the best fitting value 
for that cosmology. 



Cosmology 


Redshift bin 


r 


log(P(k=0.1)) 


EdS 


0.3 < z < 2.2 


0.16 


3.66 ±0.03 


EdS 


0.3 < ^ < 1.4 


0.16 


3.65 ±0.05 


EdS 


1.4 < z < 2.2 


0.16 


3.67 ±0.05 


A 


0.3 < z < 2.2 


0.13 


3.90 ±0.03 


A 


0.3 < z < 1.4 


0.13 


3.88 ± 0.04 


A 


1.4 < 2 < 2.2 


0.13 


3.93 ± 0.04 



at 1.4 < z < 2.2. One possible, though unlikely explanation 
for this observed flattening could be evolution of a scale- 
dependent QSO bias on ^ 100h~^Mpc scales. If this were 
the case it could be a major problem for the cosmological 
parameter analysis in the next section, and other similar 
analyses that use galaxies or QSOs as tracers of large-scale 
structure. Alternatively, the flattening could be due to the 
different shape of the window function in the two redshift 
bins. However, the two window functions are very similar, 
and if anything, the high redshift window function is slightly 
more compact, so we would expect less suppression of large- 
scale power, not more, at high redshift. 

If our assumption of an EdS cosmology r{z) were in- 
correct, as suggested by the low measured values of F and 
hence fl,nh (sec Section 5) then we would have a simple ex- 
planation for the observed flattening. If, for example, the A 
cosmology were the true cosmology, then our assumption of 
an EdS r{z) would lead to an underestimate of the comov- 
ing distance scale, by about 25 per cent in the low redshift 
bin, and 50 per cent in the high redshift bin. If we apply 
this correction to distance scale of the power spectrum de- 
termined assuming EdS, we then obtain consistent values of 
F = 0.12 ± 0.03 at 0.3 < z < 1.4 and F = 0.13 ± 0.03 at 
1.4 < z < 2.2. The observed evolution in shape assuming 
an EdS r{z) could therefore be taken as further marginal 
evidence against an EdS cosmology. 



5 COMPARISON WITH MODELS OF LARGE 
SCALE STRUCTURE 

5.1 The Hubble Volume simulation 

The Hubble Volume simulation was introduced in section 2.2. 

The power spectrum determined from biased mock QSO 
catalogues, created using the final 2QZ selection function 
imprinted on the Hubble Volume simulation, is shown in 
Fig. 10, where it is compared to the real space input mass 
power spectrum for the simulation, together with the 2QZ 
power spectrum with A assumed. The simulation was de- 
signed to mimic the 2QZ as closely as possible, and so any 
differences in shape should give a good test of the cosmolog- 
ical assumptions. 

The first thing to note is that over a wide range 
of scales, the output rcdshift-space mock catalogue power 
spectrum follows very closely the real space input mass 
power spectrum for the simulation. On the largest scales 
(log(A;//iMpc~^) ^ —1.5), this demonstrates that the effects 
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Figure 9. The power spectrum of QSOs measured at difFereiit 
redshifts. The EdS r(z) is assumed in the upper plot and A is 
assumed in the lower plot. The open cireles show the power spec- 
trum of QSOs with redshifts in the range 0.3 < z < 1.4 and the 
filled circles show the power spectrum of QSOs with redshifts in 
the range 1.4 < z < 2.2. The errors are 1<t FKP errors. The points 
are slightly offset for clcirity. 



of the 2QZ window, which were also imprinted on the Hub- 
ble Volume simulation, are relatively small, as shown in sec- 
tion 3.1. This also shows that at scales of log(fc//i Mpc^"') ^ 
— 1.0 the biasing prescription we chose to create the mock 
catalogues (see Hoyle et al. (2002) for details) has no effect 
on the shape of the power spectrum. On smaller scales, how- 
ever, we sec a slight reduction in the output power from the 
simulation, relative to the real space input power spectrum. 
This can be attributed to redshift-space distortions caused 
by small-scale peculiar velocities of the mock QSOs. 

A comparison between the Hubble Volume and 2QZ 
power spectra shows that the simulation provides a very 
good fit to the observed data over virtually the whole range 
observed. The 2QZ power spectrum is therefore entirely con- 
sistent with the ACDM cosmological model assumed to cre- 
ate the Hubble Volume. Whilst performing detailed analyses 
of N-body simulations such as the Hubble Volume is very 
useful, providing a good test of the power spectrum esti- 
mator and a way of examining possible systematic effects. 



Figure 10. The points show the 2QZ power spectrum with A 
assumed. The dashed line shows the real space input power spec- 
trum for the Hubble Volume simulation. The solid line shows the 
Hubble Volume ACDM mock QSO catalogue redshift space power 
spectrum (see Hoyle et al. 2002), updated to match the final 2QZ 
selection function , with la errors indicated by the dotted lines. 

it is also very time consuming. We need to compare many 
different CDM models to derive constraints on cosmological 
parameters. Therefore in the next section we instead com- 
pare the QSO power spectrum to model CDM power spectra 
with a variety of different cosmologies. 

5.2 Fitting model power spectra 

Assuming that QSO bias is not scale-dependent on the scales 
we are probing, the shape of the QSO power spectrum pro- 
vides new constraints on the matter and baryonic contents of 
the Universe at a different (much higher) redshift than that 
of the 2dFGRS, or other galaxy / cluster power spectrum 
determinations. Therefore, the QSO power spectrum mesr 
surement at z ~ 1.4 can provide a powerful new test for cos- 
mological models that have previously only been constrained 
locally {z ^ 0) and at recombination {z ~ 1000). Model 
power spectra with a variety of CDM cosmologies have been 
created using transfer function formulae from Eisenstein & 
Hu (1998) to compare with the 2QZ data. We have chosen 
to vary ili/Qrn (where 0„, = ficdm + ^t), and Q,mh, con- 
sidering the range < Q,b/Q,m < 0.5 and < ^Imh < 1.0 in 
order to determine constraints on these important parame- 
ters. We have assumed a scale invariant spectrum of initial 
fluctuations. 

The models generated have been numerically convolved 
with the 2QZ window function. We noted in the comparison 
with the Hubble Volume simulation that small-scale power is 
slightly suppressed in the redshift-space power spectrum by 
small-scale peculiar velocities of the QSOs. A second sim- 
ilar effect also affects the 2QZ sample. There is an uncer- 
tainty in determining QSO redshifts from the low resolution, 
low S/N 2dF spectra of 6z ~ 0.0035 (Croom et al. 2001a), 
which, at the average redshift of the survey, assuming A, 
corresponds to an uncertainty in the line-of-sight distance of 
~ 5 Mpc. This introduces an apparent velocity disper- 
sion of CTj, ~ 600 km s~^, which can be added in quadrature 
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Figure 11. The plots show the best fit CDM models (solid lines) 
to the 2QZ power spectra (points). In the upper plot, the EdS 
r{z) is assumed. The solid line shows the best fit CDM model with 
nf,/fl^ = 0.18 and Umh = 0.24. The dotted line shows a CDM 
model from the second branch of solutions (Oj,/f2rn = 0.42 and 
Clmh = 0.70) that fit the 2dFGRS power spectrum. The lower 
plot shows the best fit CDM model to the 2QZ power spectra 
assuming the A r{z), with Qd/Qm = 0.18 and Qmh = 0.19 (solid 
line). The dotted line shows, for comparison, the same model prior 
to correction for window function and redshift-space distortion 
effects. 

to the true QSO small-scale velocity dispersion. We can ap- 
ply a simple model to correct the power spectrum for these 
redshift-spaco distortions (Ballingcr ct al. 1996). Following 
Outram et al. (2001) we assume an apparent QSO line of 
sight pairwise velocity dispersion, due to both redshift deter- 
mination errors and peculiar velocities, of 700 km s~^, and a 
value of /3 w /h = 0.39 to describe the coherent peculiar 
velocities on largo, linear scales due to the infall of galaxies 
into overdense regions. We do not attempt a correction for 
geometric distortions introduced by incorrect cosmological 
model assumptions when measuring the power spectrum. 

For this analysis, we have assumed that the data points 
are independent. Duo to the compactness of our window 
function, this should be a good approximation. We have 
also had to assume a cosmology to derive the QSO r{z) 
prior to measuring the QSO power spectrum. Ideally one 
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Figure 12. Filled contours of decreasing likelihood in the Qm.h 
— Qf,/(lrn plane, marginalizing over h and the power spectrum 
amplitude, for fits to the 2QZ power spectra, assuming EdS 
for the upper plot, and A for the lower plot. Countours axe 
plotted for a one-parameter confidence of 68 per cent (dashed 
contour), and two-parameter confidence of 68, 95 and 99 per 
cent (solid contours). + marks the best fit models to the 2QZ 
data (Qi,/Qrn = 0.18, Qmh = 0.24, or with lower likelihood, 
Qlj/Qm = 0.42, Qmh = 0.70 assuming EdS, and Qb/Qm = 0.18, 
f^mh = 0.19 assuming A), and X marks the best fit models 
(Hb/f^m = 0.18, n^h = 0.23 assuming EdS, and n^/nm = 0.15, 
^mh = 0.2 assuming A) determined in a similar analysis on the 
2dFGRS data (Percival et al. 2001). 

would like to fit each power spectrum model to the power 
spectrum measured assuming the appropriate cosmology. 
However, this would be very time consuming, and beyond 
the scope of this paper. Therefore we limit our analysis 
to power spectra measured assuming either an EdS or A 
r{z). Despite assuming a particular cosmology to derive the 
r{z), we do not initially restrict the fitted parameters to 
match these cosmologies. We determine the likelihood of 
each model by calculating the value of each fit, per- 
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Figure 13. Filled contours of decreasing likelihood in the Qmh 
— Qi,/Qrn plane, marginalizing the power spectrum amplitude 
but fixing h = 0.7, for fits to the 2QZ power spectra, assuming 
EdS for the upper plot, and A for the lower plot. Contours are 
plotted for a one-parameter confidence of 68 per cent (dashed 
contour), and two-parameter confidence of 68, 95 and 99 per 
cent (solid contours). -|- marks the best fit model to the 2QZ 
data {Qi,/Qm, = 0.12, Qmh = 0.21, or with lower likelihood, 
n^/nm = 0.43, n,nh = 0.82 assuming EdS, and nt/Um = 0.17, 
^mh = 0.19 assuming A), and X marks the best fit models 
(nt/nm = 0.18, = 0.23 assuming EdS, and Clb/^m = 0.15, 

^mh = 0.2 assuming A) determined in a similar analysis on the 
2dFGRS data (Percival et al. 2001). In each case, the dotted ver- 
tical line shows the value of Qmh assumed to calculate r{z). 



formed over the range -1.9 < log(fc//iMpc-^) < -0.75, 
and assuming that C oc exp(— x^/2). For each model, we 
allowed the amplitude to vary (to account for linear QSO 
bias; we assumed that QSO bias is scale independent on 
scales log(fc//i Mpc^^) < —0.75), choosing the value that 
maximised the likelihood of the fit. 

We also allowed h to vary, within the range 0.4 < h < 
0.9, choosing the value that maximised the likelihood of the 



fit. As the shape of the power spectrum is primarily depen- 
dent on Q,mh, and relatively invariant to changes in h for a 
given iimh, the value of h has relatively little effect on the 
results. Repeating the fit, with h fixed at ft = 0.7, following, 
for example, the results of Tanvir ot al. (1995) and the HST 
Key Project (Mould et al. 2000) only marginally improved 
the constraints on the parameters of interest. 

The best fit models, shown in Fig. 11, fit the data well 
over the full range of scales considered. The effect of the win- 
dow function and redshift-space corrections can be seen in 
the A plot. The solid line in the plot shows the best fit model 
to the A 2QZ power spectrum, with Qf,/flm = 0.18 and 
flmh = 0.19. The dotted lino shows the same model prior 
to correction for window function and redshift-space distor- 
tion effects. The effect of the window function can be seen 
on large scales, slightly supressing the power, and on small 
scales the effect of redshift-space distortions can be seen. 
This lack of small scale power steepens the power spectrum, 
loading to a lower value of F. If not corrected for, therefore, 
it would lead to an under-estimate of flmh. The effect of 
the window function is slightly larger when an EdS r(z) is 
assumed. 

Likelihood contours in the Qmh - Qt/Qm plane for fits 
to the 2QZ power spectra are shown in Fig. 12 (allowing 
ft to vary), and Fig. 13 (with fixed ft = 0.7). Assuming 
an EdS QSO r{z), we find a degeneracy between the two 
parameters. The best fit model heis Qt/Qm ~ 0.18 ± 0.09 
and Clmh = 0.24 ± 0.04, however a second solution with 
Qb/Qm = 0.42 and flmh = 0.70 cannot be rejected at 95 per 
cent confidence from this analysis alone. For comparison, we 
show the best fit derived from the 2dFGRS data (Percival et 
al. 2001), which is in good agreement with the result derived 
hero. The degeneracy we sec is also very much like that found 
in the 2dFGRS data. We note that the best fit value of Qmh 
obtained is totally inconsistent with the assumed EdS cos- 
mology. Whilst the second solution is consistent with EdS, it 
is strongly rejected by other analyses such as estimates of the 
baryon content from big-bang nucleosynthesis (O'Mcara et 
al. 2001) and recent Cosmic Microwave Background (CMB) 
results (e.g. Sievers et al. 2002). 

The best fit model to the A r{z) 2QZ data has Qt/Qm = 
0.18±0.10 and Qmh = 0.19±0.05. Although the contours are 
again elongated in the direction of the degeneracy seen in the 
EkiS plot, the second peak seen in the greyscale is ruled out 
at over 95 per cent confidence. Due to the huge volume of the 
2QZ, particularly in a A cosmology, we are able probe larger 
scales than with the 2dFGRS. Whilst the second solution 
is not a bad fit on scales log(A;/ftMpc~^) > —1.5, we can 
discriminate between the two models on the largest scales 
probed by the 2QZ. When repeating the fits, considering 
only models with ft = 0.7, we see a marginal tightening 
of the contours, shown in Fig. 13, with best fit models of 
Qi/Qm = 0.17±0.09 and Qmh = 0.19±0.04. This is totally 
consistent with the A cosmology assumed to calculate r{z). 

These results agree well with those obtained by Ya- 
mamoto (2002) who produced a similar analysis of the cos- 
mological parameter constraints that could be derived from 
the 2QZ 10k catalogue. The values of Qmh and Qh/Qm ob- 
tained assuming either r{z) arc very similar, indicating that 
these results are largely independent of our r{z) assumption. 
When fitting models in the Qmh — Qb/Qm plane we have 
considered only a family of CDM models with a scale invari- 
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ant spectrum of initial fluctuations (n = 1). If we relaxed 
these assumptions we would introduce new degeneracies (Ef- 
stathiou et al. 2002), and significantly reduce the constraints 
on the above parameters. However, the degeneracies between 
cosmological parameters derived from power spectral anal- 
yses of large-scale structure (e.g. probed by the 2QZ or the 
2dFGRS) aJone are different from those from other datasets, 
such as the CMB, and Efstathiou et al. demonstrated that 
by combining such datasets, these degeneracies can be bro- 
ken, leading to strong cosmological constraints. 



6 CONCLUSIONS 

We have presented a power spectrum analysis of the fi- 
nal 2dF QSO Rcdshift Survey catalogue containing 22652 
QSOs. Utilising the huge volume probed by the QSOs, we 
have accurately measured the power over the range — 2 < 
log(fc/^t Mpc^"*^) < —0.7, thus covering over a decade in 
scale, and reaching scales of ~ 500/i^^Mpc. The A r{z) QSO 
power spectrum can be well described by a model with shape 
parameter F = 0.13 ± 0.02. If an Einstein-de Sitter model is 
instead assumed, a slightly higher value of F = 0.16 ± 0.03 
is obtained. These low measurements of F indicate signif- 
icant large-scale power and assuming either cosmology we 
can strongly rule out a turnover in the power spectrum at 
scales below ~300/i"^Mpc (log(fe//iMpc"^) ~ -1.7). The 
amplitude of clustering of the QSOs at a ~ 1.4 is similar 
to that of present day galaxies, and an order of magnitude 
lower than present day clusters. A comparison of the power 
spectra of QSOs at high and low redshift shows little evi- 
dence for any evolution in the amplitude of clustering. 

A comparison with the Hubble Volume ACDM simu- 
lation shows very good agreement over the whole range 
of scales considered. To derive constraints on the matter 
and baryonic contents of the Universe, we fit CDM model 
power spectra (assuming scale-independent bias and scale- 
invariant initial fluctuations), convolved with the survey 
window function, and corrected for rcdshift space distor- 
tions. Assuming a A cosmology r{z), we find that mod- 
els with baryon oscillations are slightly preferred, with the 
baryon fraction Qi,/n,n ~ 0.18 ± 0.10. The overall shape 
of the power spectrum provides a strong constraint on Qmh 
(where h is the Hubble parameter), with Qmh = 0.19±0.05. 
These new constraints are derived at 2 ~ 1.4, which corre- 
sponds to a look-back time of approximately two-thirds of 
the age of the Universe. Cosmological constraints at this 
intermediate epoch provide a strong test for models previ- 
ously constrained to match CMB observations at z ~ 1000 
and the local Universe at z ~ 0. 
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